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Abstract 

We discuss the status of the Gribov-Pontecorvo (GP) solution to the solar 
neutrino problem. This solution naturally appears in bimaximal neutrino 
mixing and reduces the solar and atmospheric neutrino problems to vacuum 
oscillations of three active neutrinos. The GP solution predicts an energy- 
independent suppression of the solar neutrino flux. It is disfavoured by the rate 
of the Homestake detector, but its statistical significance greatly improves, 
when the chlorine rate and the boron neutrino flux are slightly rescaled, and 
when the Super-Kamiokande neutrino spectrum is included in the analysis. 
Our results show that rescaling of the chlorine signal by only 10% is sufficient 
for the GP solution to exist, if the boron-neutrino flux is taken 10 - 20% 
lower than the SSM prediction. The regions allowed for the GP solution in 
the parameter space are found and observational signatures of this solution 
are discussed. 



I. INTRODUCTION 



Vacuum oscillations of maximally mixed and v T neutrinos with (Am 2 ) atm ~ 3 x 
1CT 3 eV 2 are the favourite explanation of the atmospheric neutrino anomaly. A natural gen- 
eralization is bimaximal mixing 
by the following matrix: 



TUf of three active neutrinos, when mixing is described 
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In the case of the mixing matrix given by Eq. (U), the solar neutrino oscillation is also 
maximal. To see this, it should be noted from Eq. (|1|) that 
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while the other orthogonal combination of these states can be considered as a new field 
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Using the above equations, one obtains 
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From Eq. (^) it follows that v e and v' = l/v^2(^ At + v T ) are a maximally mixed pair, and the 
flavour eigenstate v e is oscillating on the way from the Sun into z/, the coherent mixture of 
and v T . 

The above exercise is relevant to the Gribov-Pontecorvo (GP) Jll| solution of the solar 
neutrino problem combined with atmospheric Vn—v T oscillations. Following ||11||, the defini- 



tion of the GP solution can be given by two conditions: (i) smallness of the oscillation length 
l v with respect to the mean distance between the Sun and the Earth (r) = L 
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< L = 1.5 x 10 
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cm, 



(5) 



and (ii) smallness of the matter corrections (MSW |T!|) in the Sun and in the Earth P0| (in 
Ref. |]TTJ] only vacuum oscillation is considered). 

Indeed, in this case the averaged survival probability for v e is 



(P ee ) = l-sin 2 2fl(sin 2 ^> 



1 - - sin 2 29 = P G p, 



(6) 



and from comparison with experimental data, P ee ~ 0.5, we come to 9 ~ 7r/4, or to bimaxi- 



mal mixing if za, 



explains atmospheric neutrino oscillations. 



Three remarks are immediately in order. 

(i) There is no theoretical reason for bimaximal mixing to be exact, and more generally 
one should consider near- bimaximal mixing P|lr5[|. 

(ii) The smallness of the matter correction effects, which we included in the definition of 
the GP solution is actually not needed in the case of near-bimaximal mixing if only matter 
effects in the Sun are included. For the exact maximal mixing, v e v' conversion in the Sun 
does not change the total survival probability (P ee ) at the surface of the Earth [15|fl. In the 
arbitrary case where the mixing between u e and v' is described by the mixing angle 9, and 
the MSW effect in the Sun converts v e into v s = cos 6sV\ + sin^z^ (subscript S here refers 
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to the surface of the Sun), the total survival probability (P ee ) = K^l^e)! 2 on the surface of 
the Earth can be readily calculated as 

(Pee) = cos 2 9 S cos 2 9 + sin 2 9 S sin 2 9 (7) 

For exact bimaximal mixing cos 2 6 = sin 2 9 = 1/2, survival probability (P ee ) = 1/2 and 
thus it does not depend on 9s, i-6. on how v e is converted in the Sun. For near-bimaximal 
mixing, Eq. (^) determines a narrow range of 9 near 7r/4, where (P ee ) is practically energy 
independent, i.e. it does not depend on 9s- The matter effect in the Earth, however, changes 
this conclusion as we shall see in the next section. 

(in) The observational data (see Fig. [I]) do not support (P ee ) being exactly energy- 
independent. While the recent Super-Kamiokande (SK) data agree well with (P ee ) being 
an energy- independent constant in the energy interval 5-14 MeV, the values (P ee ) from 
three different experiments, GALLEX-GNO/SAGE, Homestake and Kamiokande/SK, are 
not exactly the same. 

The aim of this paper is to discuss quantitatively the status of the GP solution. An 
interesting region in the parameter space is given by Am 2 < 10~ 3 eV 2 and 9 ~ 7r/4, where 
bimaximal mixing is characterized by (Am 2 ) atm and (Am 2 ) so i not very different from each 
other, and where the MSW effect is allowed as a small correction. 

Oscillations with energy-independent suppression were suggested and studied in many 



works before |i2|-4T6| and most notably in the recent work of Ref. jD]]. In all these works 



the authors have realized that the observed rate in the chlorine experiment (Homestake 



25fl ) contradicts the energy-independent suppression, and it has to be taken larger than the 
observed one for it to work. In this paper we argue that 10% excess could be sufficient. Note, 
that the solution with energy-independent suppression is more general than the GP solution, 
because in the latter l u <C L is assumed, while the energy- (quasi) independent solution might 
appear in some other regions of the parameter space. 



II. PARAMETER SPACE REGIONS FOR THE GP SOLUTION 

In this section we shall calculate the regions allowed for the GP solution in the parameter 
space Am 2 , tan 2 9. We first define the oscillation parameter space where the solar neutrino 
survival probability behaves effectively as the GP one. In order to do so we impose the 
condition that for any of the i solar neutrino fluxes (integrated over the different production 
point distributions) the survival probability in the relevant range of energies does not differ 
by more than 10% (1%) from P GP given by Eq. (^|): 

\PUE/Arn*,9)-P G Ae)\ < Q1 {QM) for >E> ^ (g) 

PsplyJ 

where E itmax and E itmin determine the range of energies in which the flux i is detected in 
present experiments. For instance, for i = pp, E PPtmin = 0.233 MeV and E PPtmax = 0.42 MeV. 
In the evaluation of the corresponding survival probabilities, we have included the matter 
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effects when propagating in the Sun and in the Earth as well as the distance interference 
term: 



PL = + P e^e + 2^1^1^008 , (9) 

where P' t A is the probability for the v e to exit the Sun in the mass eigenstate i, while Pf e is 
the probability for the mass eigenstate % arriving at the Earth to reach the detector as a v e . 
L is the average distance between the Sun and the Earth. 

In Fig. ^ we show the parameter space Am 2 , tan 2 9 where the condition given by Eq. (§) 
is verified at 10% (lighter shadow) and 1% (darker shadow). The only interesting sector of 
the effective-GP region in this parameter space is located at large Am 2 around the maximal 
mixing line tan 2 9 = 1, where matter effects in the Sun are suppressed. This region is limited 
from above by the CHOOZ reactor data f27fl. Only in this sector there is an overlap with 
the rate- and spectra-allowed regions (see below). 

As was discussed above, for maximal mixing the matter effects in the Sun do not alter the 
energy-independent survival probability P ee on the way from the production point inside the 
Sun to the surface of the Earth. However Earth matter effects make P ee energy- dependent 
in the regions of maximal mixing at 10~ 5 iS Am 2 < 10~ 8 . In contrast with our calculations, 
the region Am 2 < 1CT 7 is found in ref. [17] as energy-independent one. We explain this 
discrepancy by two effects: 

• For 1CT 8 < Am 2 < 10 -7 Earth matter effects for pp-neutrinos result into an energy 
dependence of the survival probability beyond 10%. 

• At Am 2 < 10~ 8 the L dependent interference term in Eq. (jfj) gives strong energy 
dependence of the 8 B flux. This term was not included in the calculations of Ref. [|17] . 



As we mentioned above, the GP solution is incompatible with the central value of the 
rate measured by the Homestake detector Rqi = 2.56 SNU fl25| . Following the prescription 



of many works, we shall use the rescaled rate R^f = 2.56/ci SNU, assuming fa > 1 to 
be a free parameter. In Fig. § we plot the y 2 function from the analysis of the three 
observed rates as a function of the fc\ factor for different constant values of the survival 
probability P GP = 0.5,0.59 and 0.71. The upper left panel corresponds to oscillations into 
active neutrinos while the lower one into sterile neutrinos. The differences between these 
two scenarios arise from the absence of NC contribution to the SK rate in case of oscillations 
into sterile neutrinos. From this figure we see that the best GP-like solution corresponds to 
survival probability slightly larger than that for maximal-mixing case (close to 0.59) both 
for active and sterile neutrinos. The quality of these solutions are considerably improved 
when allowing a 30-50 % increase in fc\- This improvement is more significant for the case 
of sterile neutrinos since the corresponding survival probability at SK agrees better with the 
data from gallium detectors (see Fig.|TJ) . 

This behaviour is also illustrated in Fig. [| where we show the Am 2 , tan 2 9 regions allowed 
by the statistical analysis of the rates of GALLEX-GNO/SAGE [H,||,|l|], SK H and 



Chlorine |25[] experiments for different values of /ci in case of active and sterile neutrinos, 
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and for Bahcall-Pinsonneault (BPOO) |26| fluxes. The solutions, following the standard 
statistical analysis (for details see Ref. [Pf| ), are shown at 99% CL. The effective GP solutions 
are marked as dark areas. Notice that they appear at fci ~ 1.1 (1-2) for active (sterile) 
oscillations and that all regions displayed have a cut at Am 2 m 8 x 10~ 4 as a consequence 
of the CHOOZ |7| bound. 

Inclusion of SK data on the energy spectra of boron neutrinos improves the quality of 
the GP solution. In Fig. |5| we display, for different values of fci, the regions allowed by 
the analysis of the rates and day-night spectrum of boron neutrinos measured by SK . 
Again the solutions are shown at 99% CL and the effective GP solutions are marked as dark 
areas. From these figures one can see that the inclusion of the spectra data results in the 
appearance of allowed regions for the GP solutions at smaller values of fci- This is a natural 
result, because the rates of GALLEX-GNO/SAGE and Homestake can be also considered as 
information about the solar neutrino spectrum, in its low energy part, however in contrast to 
the low energy part of the spectrum, the GP solution describes well the spectrum observed 
in SK. 

In Figs. |]-|5] we have used the boron neutrino flux as calculated in the Standard Solar 
Model PE| ]. This flux has a large theoretical uncertainty mostly due to uncertainties in 
the pBe cross-section. In order to study the effect of a possible deviation of the 8 B flux 
from the SSM prediction |]26|| , we shall introduce the rescaled boron neutrino flux defined as 
$b = /b$| SM with $| SM given by |26]. In the central panels in Fig. |3], we plot the y 2 function 
from the analysis of the three observed rates as a function of fci for different constant values 
of the survival probability Pgp when the factor /b is left free. The upper central panel 
corresponds to oscillations into active neutrinos, the lower one into sterile neutrinos. In the 
right panels we show the corresponding values of /b, which give the best agreement with 
the data for each value of fci and Pgp- From this figure we find that although using a free 
/b < 1 leads to a further increase of the statistical significance of the GP solution, it has a 
smaller impact than the corresponding variation of /ci- It, however, allows for the presence 
of GP solutions with smaller fci- This is particularly the case for oscillations into active 
neutrinos. 

In Fig. ^we plot the allowed regions from the analysis of the rates and day-night spectrum 
of 8 B neutrinos measured by SK for different values of fci and /b. For the sake of concreteness 
we have chosen the /b factor that gives a better fit to the three rates for each value of fci for 
maximal mixing. Namely, /b is chosen as /b ~ 0.8 (0.9) for oscillations into active (sterile) 
neutrinos. In Fig. |6| the left panels correspond to oscillations into active neutrinos and the 
right ones into sterile neutrinos. Comparing this figure with the corresponding panels in 
Figs. |5| we see that lowering the 8 B normalization leads to a larger overlap between the 
allowed LMA region and the GP solution already for fci < 1.1, in the case of oscillations 
into active neutrinos. 



III. CONCLUSIONS 

It could be that nature has chosen the most unsophisticated scheme of neutrino oscilla- 
tions: three active neutrinos with (nearly) bimaximal mixing. Mixing of and v T explains 
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the atmospheric neutrino anomaly, and of v e and v' = (l/v / 2)(^ + v T ) the solar neutrino 
deficit. In this case, the GP solution, provided by condition (|5|), naturally appears, and it is 
characterized by an energy- independent survival probability (P ee ) — Pgp- 

The GP solution is disfavoured by the Homestake rate, but describes well the other rates 
as well as the energy spectrum observed in SK. The statistical significance of the GP solution 
strongly improves if one assumes rescaling of the chlorine rate by a factor fc\ = 1.1—1.5, while 
some further improvement arises if the 8 B neutrino flux is also rescaled by a factor /b = 0.7- 
0.9. In particular if the 8 B flux happens to be 10-20% lower than the BPOO-predicted central 
value, the GP solution for maximal mixing in active oscillations would be a llowed with a 
chlorine rescaling factor fci ^5 1.1. 

The GP solution will be directly searched for in the KamLand experiment . Detection 
of reactor neutrinos can result in the measurement of Am 2 in the interval 10~ 3 -3 x 10 -6 eV 2 
for large mixing ang les. If Am 2 is found outside the LMA MSW region or inside it at 
9 ~ 7r/4, would mean the discovery of the GP solution. 

In low energy solar neutrino experiments the signatures of the GP solution are the or- 
dinary suppression of 7 Be neutrinos given by a factor Pqp = \ sin 2 29 and the absence of 
anomalous seasonal variations (beyond the geometrical ones). These features can be clearly 



seen in Borexino [|30| and KamLand p9| experiments. 
Note 

This work was presented by M.C.Gonzalez-Garcia at the Gran Sasso Laboratory at the 5th 
Topical Workshop on "Solar Neutrinos: Where are the Oscillations?" (March 2001). On 
March 29 the preprint by S. Choubey, S. Coswami, N. Gupta and D.P. Roy |L7j appeared 
in the net. The basic assumptions they used, rescaling of the chlorine and boron fluxes, are 
the same as in our paper, but we are considering a GP solution that is not, in principle, 
identical to the energy-independent solution, studied in the aforementioned paper. The most 



noticeable difference in our results relates to the low Am 2 solutions found in Ref. Wfi] and 



shown in their figures 1 -4. They are not present in our solutions partly due to interference 



term given in our Eq.(pf) and disregarded in Ref. |IT| 
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FIG. 1. Ratios of observed rates to the BPOO prediction for the existing experiments. In the 
case of oscillations the ratios for all experiments, except Super-Kamiokande, are equal to the v e 
survival probability P ee . For Super-Kamiokande the ratio SK gives P ee in the case of oscillation 
to sterile neutrinos, while for the case of oscillation to active neutrino we plot the corresponding 
p 

1 ee,ac- 




tan 2 i? 

FIG. 2. Regions in Am 2 , tan 2 6 parameter space, where the v e survival probability P ee (E) 
differs from the energy-independent GP survival probability Pgp by less than 10% (lighter shadow) 
and less than 1% (dark shadow). See the text for details. 
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FIG. 3. x 2 analysis of the three observed rates as a function of the fci factor for different 
constant values of the survival probability Pgp = 0.5,0.59 and 0.71 (solid, dashed and dotted lines 
respectively). The upper (lower) panels correspond to oscillations into active (sterile) neutrinos. 
In the left panels we have used BP00 boron fluxes (/b = 1), while in the central panels /b is left 
free to optimize x 2 . The resulting /b values which optimize x 2 f° r given fci and Pqp are plotted 
in the right panels. 
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FIG. 4. 99% CL regions allowed by the analysis of the experimental rates in 
GALLEX-GNO/SAGE, SK and Chlorine experiments, for BPOO fluxes, for different values of 
/ci and for oscillations into active (left panels) and sterile (right panels) neutrinos. The dots mark 
the position of the best-fit points in each panel. The effective GP solutions are marked as dark 
areas. 
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FIG. 5. 99% CL regions allowed by the analysis of the three experimental rates and the SK 
day-night spectrum, for BPOO fluxes, for different values of fci and for oscillations into active (left 
panels) and sterile (right panels) neutrinos. The dots mark the position of the best-fit points in 
each panel. The effective GP solutions are marked as dark areas. 
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FIG. 6. 99% CL regions allowed by the analysis of the three experimental rates and the 
SK day-night spectra for different values of fci and /b, for oscillations into active neutrinos (left 
panels) and sterile neutrinos (right panels). The effective GP solutions are marked as dark areas. 
The dots mark the position of the best-fit points in each panel. 
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